Kinetic and thermodynamic analysis of proteinlike heteropolymers: Monte Carlo histogram technique

Using Monte Carlo dynamics and the Monte Carlo Histogram Method, the simple three-dimensional 27 monomer lattice copolymer is examined in depth. The thermodynamic properties of various sequences are examined contrasting the behavior of good and poor folding sequences. The good (fast folding) sequences have sharp well-defined thermodynamic transitions while the slow folding sequences have broad ones. We find two independent transitions: a collapse transition to compact states and a folding transition from compact states to the native state. The collapse transition is second order-like, while folding is first order. The system is also studied as a function of the energy parameters. In particular, as the average energetic drive toward compactness is reduced, the two transitions approach each other. At zero average drive, collapse and folding occur almost simultaneously; i.e., the chain collapses directly into the native state. At a specific value of this energy drive the folding temperature falls below the glass point, indicating that the chain is now trapped in local minimum. By varying one parameter in this simple model, we obtain a diverse array of behaviors which may be useful in understanding the different folding properties of various proteins.Comment: LaTeX, 16 pages, figures in separate uufile. Requires psfig.sty Minor revision, fixed typo in preprint number (no other changes

Folding Kinetics of Protein Like Heteropolymers

Using a simple three-dimensional lattice copolymer model and Monte Carlo dynamics, we study the collapse and folding of protein-like heteropolymers. The polymers are 27 monomers long and consist of two monomer types. Although these chains are too long for exhaustive enumeration of all conformations, it is possible to enumerate all the maximally compact conformations, which are 3x3x3 cubes. This allows us to select sequences that have a unique global minimum. We then explore the kinetics of collapse and folding and examine what features determine the various rates. The folding time has a plateau over a broad range of temperatures and diverges at both high and low temperatures. The folding time depends on sequence and is related to the amount of energetic frustration in the native state. The collapse times of the chains are sequence independent and are a few orders of magnitude faster than the folding times, indicating a two-phase folding process. Below a certain temperature the chains exhibit glass-like behavior, characterized by a slowing down of time scales and loss of self-averaging behavior. We explicitly define the glass transition temperature (Tg), and by comparing it to the folding temperature (Tf), we find two classes of sequences: good folders with Tf > Tg and non-folders with Tf < Tg.Comment: 23 pages (plus 10 figures included in a seperate file) LaTeX, no local report nu

Genome-wide mapping of the distribution of CarD, RNAP σA, and RNAP β on the Mycobacterium smegmatis chromosome using chromatin immunoprecipitation sequencing

CarD is an essential mycobacterial protein that binds the RNA polymerase (RNAP) and affects the transcriptional profile of Mycobacterium smegmatis and Mycobacterium tuberculosis [6]. We predicted that CarD was directly regulating RNAP function but our prior experiments had not determined at what stage of transcription CarD was functioning and at which genes CarD interacted with the RNAP. To begin to address these open questions, we performed chromatin immunoprecipitation sequencing (ChIP-seq) to survey the distribution of CarD throughout the M. smegmatis chromosome. The distribution of RNAP subunits β and σA were also profiled. We expected that RNAP β would be present throughout transcribed regions and RNAP σA would be predominantly enriched at promoters based on work in Escherichia coli [3], however this had yet to be determined in mycobacteria. The ChIP-seq analyses revealed that CarD was never present on the genome in the absence of RNAP, was primarily associated with promoter regions, and was highly correlated with the distribution of RNAP σA. The colocalization of σA and CarD led us to propose that in vivo, CarD associates with RNAP initiation complexes at most promoters and is therefore a global regulator of transcription initiation. Here we describe in detail the data from the ChIP-seq experiments associated with the study published by Srivastava and colleagues in the Proceedings of the National Academy of Science in 2013 [5] as well as discuss the findings from this dataset in relation to both CarD and mycobacterial transcription as a whole. The ChIP-seq data have been deposited in the Gene Expression Omnibus (GEO) database, www.ncbi.nlm.nih.gov/geo (accession no. GSE48164)

Non-equilibrium raft-like membrane domains under continuous recycling

We present a model for the kinetics of spontaneous membrane domain (raft) assembly that includes the effect of membrane recycling ubiquitous in living cells. We show that the domains have a broad power-law distribution with an average radius that scales with the 1/4 power of the domain lifetime when the line tension at the domain edges is large. For biologically reasonable recycling and diffusion rates the average domain radius is in the tens of nm range, consistent with observations. This represents one possible link between signaling (involving rafts) and traffic (recycling) in cells. Finally, we present evidence that suggests that the average raft size may be the same for all scale-free recycling schemes.Comment: 8 pages, 5 figure

Comprehensive evaluation of differential gene expression analysis methods for RNA-seq data

A large number of computational methods have been developed for analyzing differential gene expression in RNA-seq data. We describe a comprehensive evaluation of common methods using the SEQC benchmark dataset and ENCODE data. We consider a number of key features, including normalization, accuracy of differential expression detection and differential expression analysis when one condition has no detectable expression. We find significant differences among the methods, but note that array-based methods adapted to RNA-seq data perform comparably to methods designed for RNA-seq. Our results demonstrate that increasing the number of replicate samples significantly improves detection power over increased sequencing depth

MicroRNA profiling of the murine hematopoietic system

BACKGROUND: MicroRNAs (miRNAs) are a class of recently discovered noncoding RNA genes that post-transcriptionally regulate gene expression. It is becoming clear that miRNAs play an important role in the regulation of gene expression during development. However, in mammals, expression data are principally based on whole tissue analysis and are still very incomplete. RESULTS: We used oligonucleotide arrays to analyze miRNA expression in the murine hematopoietic system. Complementary oligonucleotides capable of hybridizing to 181 miRNAs were immobilized on a membrane and probed with radiolabeled RNA derived from low molecular weight fractions of total RNA from several different hematopoietic and neuronal cells. This method allowed us to analyze cell type-specific patterns of miRNA expression and to identify miRNAs that might be important for cell lineage specification and/or cell effector functions. CONCLUSION: This is the first report of systematic miRNA gene profiling in cells of the hematopoietic system. As expected, miRNA expression patterns were very different between hematopoietic and non-hematopoietic cells, with further subtle differences observed within the hematopoietic group. Interestingly, the most pronounced similarities were observed among fully differentiated effector cells (Th1 and Th2 lymphocytes and mast cells) and precursors at comparable stages of differentiation (double negative thymocytes and pro-B cells), suggesting that in addition to regulating the process of commitment to particular cellular lineages, miRNAs might have an important general role in the mechanism of cell differentiation and maintenance of cell identity

Derivation of Multipotent Mesenchymal Precursors from Human Embryonic Stem Cells

BACKGROUND: Human embryonic stem cells provide access to the earliest stages of human development and may serve as a source of specialized cells for regenerative medicine. Thus, it becomes crucial to develop protocols for the directed differentiation of embryonic stem cells into tissue-restricted precursors. METHODS AND FINDINGS: Here, we present culture conditions for the derivation of unlimited numbers of pure mesenchymal precursors from human embryonic stem cells and demonstrate multilineage differentiation into fat, cartilage, bone, and skeletal muscle cells. CONCLUSION: Our findings will help to elucidate the mechanism of mesoderm specification during embryonic stem cell differentiation and provide a platform to efficiently generate specialized human mesenchymal cell types for future clinical applications